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Abstract 

We investigate effects of a nonadiabatic electron-phonon(boson) interaction 
on the quasiparticle seff-energy in the lowest order in the coupling constant. 
Existing approaches either overestimate, or underestimate these effects be- 
cause of different approximations for momentum and frequency dependences 
of the vertex corrections. The connection between the nonadiabaticity and a 
possible instability of the interacting Fermi system is discussed as well. 

I. INTRODUCTION 

The standard description of electron-phonon interaction (EPI) in metals is based on the 
so-called Born-Oppenheimer0( or adiabatic) theorem. Its main statement says that electrons 
are not sensitive to the motion of ions and are influenced only by their static electric field. 
One can to say that "fast" electrons, due to electroneutrality, follow "slow" ions. Narrow 
bands and a strong electron-phonon interaction make their velocities comparable in order 
of magnitude. An appropriate mathematical description is the standard Feynman-Dyson 
perturbation theory, where phonons are considered to be uneffected by the EPI. However, 
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there are some processes which in higher orders violate this approximation. Mathematically 
this can be taken into account by virtue of the so-called "vertex corrections". These effects 
were considered by Migdali who showed that their contributions are small, of the order of a 
small parameter Qph/W, where is a characteric energy, defined by the bandwidth W or 
the Fermi energy Ep. 

It was noted by Migdal himsel^l that at g = the vertex corrections are not small (so they 
are important for optical conductivity and Raman scattering). This case was considered in 
detail by Engelsberg and Schriefferi. Migdal's theorem is also violated for a one-dimensional 
Fermi surfacei. 

Traditionally nonadiabatic corrections were considered to be small because the Migdal 
parameter Qph/W was smalliffli'i. But recently some materials (for example, fullerenes 
and high-Tc superconductors) where Qph ~ W were discovered. So Migdal's theorem about 
the smallness of the nonadiabatic correctiondli'i can be violated. This gives grounds for not 
only discussing these corrections, but even examining an antiadiabatic limit ilp^ ^ lyi. But 
even with small Qph/W there are contradictory conclusions about the importance of such 
corrections in the normal and superconducting state (the authors of Refs.i&i consider them 
to be negligible, while the authors o^lii^iii^ii think the opposite is true). The most intriguing 
question is whether the violation of adiabaticity leads to an increase of the interelectron 
interaction (and, as a result, to an increase of T^). There are contradictory results about 
the effect of the vertex corrections on T^. For example, Takadai reported some increase of 
Tc for small parameters Vtph/Ep. Even more significant increase was claimed in0. At the 
same time, in Refs.izlii'iS a decrease of Tr was obtained due to the vertex corrections. 



One can distinguish between the two sorts of the nonadiabatic corrections to the electron 
self-energy: 1) the vertex ones, resulted from the higher order corrections of perturbation 
theory, as vertex function F 7^ 1 even at an infinite electron (or hole) band — 00 < e < 
+ooi'iil, and 2) the corrections due to the finite band width, -W < e < iyillHli§i. 

Nonadiabatic corrections in the superconducting state enter the expression for in two 
ways: indirectly through renormalization of Z and directly in the equation for the order 



parameter A. At the same time, analytical calculation of anomalous (crossing) diagrams 
is rather difficult. Therefore before proceeding to calculations in the superconducting state 
we want to clarify the role of such corrections in the normal state. The main goal of 
this paper is to compare different approaches to calculation of nonadiabatic corrections to 
the electron quasiparticle self-energy in the normal state (or more precisely to the mass 
renormalization factor Z). There are two methods to take nonadiabatic corrections into 
account. One of them, which can be called Migdal's one, is based on the solution to the 
Bethe-Salpeter equation for the vertex function&ifl in the ladder approximation. In the 
lowest approximation, the first correction to the unity vertex is determined by the diagram 
in Fig. 1. Then with the obtained vertex function the self-energy E is calculatedS'ffllii. This 
traditional method does not allow to take into account higher-order diagrams analytically 
because of complexity of the integration over momenta. 

In 1989 the authors of the Ref.i3 considered a non-ladder approximation, using the Ward 
identity. Later Y. Takada used the same idea in his methodi, which he called the gauge- 
invariant self-consistent (CISC) method. In this method, based on the Ward identityi 



the vector term is neglected and the scalar vertex function is chosen to be a functional of the 
self-energy, which is supposed to be independent from momentum. Recently there appeared 
another methodil, aiming at improving the CISC approach as it takes into account the 
momentum dependence of the vertex function. So, when comparing these methods with 
the Migdal's one, we also discuss the validity of these assumptions. In fact, for the self- 
energy we calculate all the diagrams up to, and including the second order. The parameter 
^ph/W = mo is taken to be small. 

For simplicity, we consider the model of the Fermi liquid with the usual electron-phonon 



tion a^FE^^) = ^XQphS{Q — 0,ph), where Qph is the phonon frequency independent from 



iojyViiuJn, iuJn — i^^f, k,k — q) — qV^iUn, ioJn — i^u, k,k — q) 
G~^{iuJn, k) - G~^{iuJn - iuj^, k - q), 



(1) 




and Einstein ph 



.ononBa spectrum with the so-called Eliashberg spectral func- 
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momentum, and A is the electron-phonon interaction constant. It is assumed that the den- 
sity of states is constant A^(0) for —W < e < W where 2W is the bandwidth. The chemical 
potential is zero, which corresponds to a half-filled band. It is assumed that A ~ 1. 

We are interested in the quasiparticle self-energy S on the imaginary axis which is deter- 
mined by the diagram in Fig.l. The phonon Green's function is DQ^iu^) = —Qpf^/{ujl + flpf^), 
where = 2z/7rT, and the electron Green's function is G{iu!n,p) = {iojn — £p — T,{iu!n,p))~^ , 
where ujn = (2n + l)7rT, is the bare electron spectrum. 

The vertex function itself (two-point one) is only of academic interest. We are, how- 
ever, interested in the quasiparticle self-energy S which affects many physical observables. 
Assuming that S(ic(j„, k) weakly depends on k, we take \k\ = Pf = const, and then 

2n 1 



S(icu„) = f dna^FE{n)Tj2 

■J 1.1 



- uji + iv(o) 

X ^ G{iu}ri - i^u, k - q)r{iun, iuJn - iuJu, q)- (2) 

The methods differ in the choice of the vertex function F, and this will be discussed 
later. One should note that r(g) enters the integral for S, and all g's contribute to it. For 
simplicity we will obtain results for T = and the variables will be changed as follows: 
z/, remaining on the imaginary axis. 



II. MIGDAL'S METHOD 

In the method which can be called Migdal's, the first correction to the unity vertex 
function F*^^^ is considered. This is shown in Figl. It has been estimated in many papersi'SQ, 
but now it is becoming especially significantEHllbecause the cases were found where as the 
parameter XVtph/W may be not small. 

The diagram in Fig.l corresponds to the expression 



^'V q' 

xG^^\iuJn — i^v')V — (f)G^'^\iuJn — i^u' — i^v,V — — O)- (3) 
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Assuming that |p — (f\ ~ pp and expanding e{j> — if — q) ^ ^{p — <f) — qVpcosO we get 
for T = 

1 + im + 1 — im 1 + im — im' 



xiln ^ + ln- 

^ l-im + — 1 



+ ln 



1 — im -\ 1 + im 1 — im + im' 

' mo ' 

1 — im + im' — Qn + — 

+ ln — (4) 

1 + im — im' + Qn + — ' 

^ ' mo 

where m = uj/Qph, m' = i^/Qph, Q = qVp/flph, = Qph/W. 

This is a nonanalytical function at uj 0, q ^ 0, i.e. the resuh depends on the order of 
taking the hmits. Let us define the dynamical and static vertex functions 

Yd = T{q = 0, z/ -> 0,uj); 

T, = Tiq-^0,u = 0,ij). (5) 

From one can get 

1 ri-^jn -I 



Tf =TM{q = 0,iy^O,uj) = X[ 



W \2 j_ ( <^ ^2J' 



1 + -^ 

rf = rM(g - o,z. = 0,^) = -A 

It is evident that the nonanalyticity mentioned above gives 

1 



-pM _ -pM _ \ 

d. ^ s ~ ^ 1 I (JiJ_\2 ■ 



At Vplfl < |z/| 



r^'^\iuj, iuj — iu, q) 

, Vtj,h UJ uj — u LO 

A— ^(arctan — arctan — arctan ■ 



V Qph Qph W + Qph 



UJ — u 



+ arctan ^ " ) - ^^^H^^) - ^(^H^^ - .g. 
-haiCLau-r— -- i — , luj 

W + ^lph IV 
which satisfies the Ward identity ([^). 
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At VfI^I ^ |z/| one has T^^\iuj,iu — iv,q) XVLph/W, which comphes with Refs.i'i. 
To calculate the self-energy we expand the last term of in Q, leaving the first two 
terms as they are. Then 

r*-^-* (icu, iu — iiy,q) = X—rr arctan 

, ijJ ijJ — V ijj ijj — V 

X (arctan — arctan — arctan — h arctan 



f22 



In Fig. 2a we show r^^)(z/) (numerical (^ and approximate analytical (^ results) for the 
case = 0, as well as the result of the numerical integration of Eq.(4). One can see that 
there is a good agreement between these plots, so the expansion in small Q is well justified. 

The self-energy (^) consists then of three terms: 

S = S« + S(2) + S(2). (8) 

The first order term isi 

Y}^^^{iijj) = —iXQ,ph arctan — h iXVtph arctan ^ — . (9) 

The vertex correction E^^^ is obtained when one substitutes G = and T = F*^^) (^ 
into (@) 

S(2)(M = Sf'^(.c.) + S(2)''(..;), (10) 

Sl2)'^(z^) =zA2-^<(m,mo,m'o), 
PfVf 

m = uj/^ph, mo = Vlph/W, m'^ = VLph/^ppVp)- 

1 7 dy 1 



(i"(m, mo, mn) = / -—-^ — arctan 

ZTT J V + 1 



271 J y"^ + I mo{m — y) 



— oo 



m m — y 
X (arctan m — arctan(m — y) — arctan — h arctan — ) 

mo ~^ rriQ ~^ 

2 1 4 

X (2 arctan ^ - -ym'o ln(—- ^ + 1)). (11) 

yrriQ 2 1/ 
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Vf vf 



where 



1 /'+°° 1 1 

dl{m, mo, rn'o) = -—ml{l + mo) / dm' arctan — — - 

2tt j-oo m'^ + 1 molm' — m) 

r 2 2m' 2 m'2 , 9 

x{ — 2 -arctan— — + In 1 + ( -f]}. (12) 

m'g mo mmQ 2 m'm^ 

For cj/fip^j < 1, fip/i/iy < 1 it gives 

4.M...=^.[^,,i;£i£,.]. (13) 

The plot of S(2) (ic(j) is shown in Fig. 2b. The numerical calculations (the exact integration 
with Eq.(^)) give similar results for the small a;'s . The existing discrepancy for larger 
frequencies is due to our approximation for F. 

Another term is the so-called rainbow diagram. As the inner part of the rainbow diagram 
corresponds to Tj^^\iuj — iu), the second order rainbow diagram equals 

2n 1 





^ Yl[G{iuJn - iuJu, k - q)]'^T.^^\iuJn - iuj^), (14) 



which results in 



(m^mo), where 



dr{m, mo) 



mo(m4 + 2m2 + 1 + ^ + 2(^)2 - ^) 

*J ^ nig ^ mo ' nig ' 



m'^ + 2'^ m ^ \/ 2 1 \ 
xfmln — ; + (arctan arctan )(m — IH ^) 

-2mln^^). (15) 

mo + 2^ ^ ' 



At uo/Vtph < 1, Qph/W < 1, it equals 



^f\ru)=r\x'u^^. (16) 
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For uo/Vtph <ti 1, ^ph/W <ti 1, ^ph/iPFVF) <^ 1 one can get 



1 + A - — A"— ^ - -A^i:^ - 



luj 8 pfVf 2 W ' 

One can see that the last expression has three terms with different denominators. When, for 



example, ppVp ~ ^ph-, the main contribution comes from the term —^X^VLph/ppVp (even 
for an infinite bandwidth W). In Fig. 2b there are shown for comparison analytical and 
numerical plots for the self energy correction S^^^. It is seen that the approximation used 
to calculate the correction is well justified. 
If one assumes W = ppVp, 

Z,,^l + \~ 2.7A^^. (17) 

One can see that these corrections lower the renormalization function. Had this tendency 
persist in all higher orders, this would have resulted in an instability (see below). 

III. CAI, LEI, AND XIE'S APPROXIMATION 

Following the paperil, the authors of Ref.il neglected the momentum dependence of the 
vertex function, considering q ^ ppu/W, and interpolated it between zero-frequency and 
infinite-frequency limits: 



\j^ + {ujn-u;n'r °o;^-|-o;2 + 2(Aoo/Ao)- 



Where Ao(^) = ^^f^, and A^(.;)/Ao(u;) = '-f'-^^ 
This results in an effective interaction 



iV(0) - ^')2 + (]2 1+ Ag^^^ 



0.293n 



Combined with the contribution from the rainbow diagram, which is the same as in the 
Migdal method (H), this gives 
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■^CLx « 1 + A - 1.38A2^. 



W' 



The non-ladder approximation which was also discussed i n0 was developed in detail by 
Takada, and this is considered below. 

IV. KOSTUR AND MITROVIC'S APPROXIMATION 

There has been a series of papers, e.g. Ref.ii, which considered the effect of the vertex 
corrections on not only for an isotropic EPI, but also for spin (antiferromagnetic) fluc- 
tuations with a pronounced scattering at the wave vector Q* = (vr,7r). In the paper^l the 
authors used the following equations for the self-energy S(/c,iu;„) = iUnW — Zn{k)]: 

rp2 



iu^nZ^^\k) = E Kn - n')X{n - n") 



n'n' 

(iek' / (iek'" / (iek"M(k, k'"; k', k")^^ 

~W J-W J-W Zn' 



X 



/// 1 / 1 //\ lUJn't'^n"'l^n"' 



Zn"Zn" 



1 

[iUjI, - {t]^, / Zn'Y][iujl„ - {ei,n / Zn»Y][iujl,„ - (ek'"/Zn"')^]' 

where \{m) = dQa'^F{Q) ^2^]^f^2 , and M(k, k'"; k', k") is a geometrical factor being a 
complicated function of momenta which is considered in detail in Ref.0. For an isotropic 
EPI and a three-dimensional spherical Fermi-surface it is approximated by its value on the 
Fermi surface M(k, k'"; k', k") = 1. It gives 

Z„ = 1 + - — - A(n - n )rKM{n, n')s„s„/a„/, where 

A(m) = A 2 la^ , cin = 2/7r arctan(iy/Z„|ci;„|), s„ = sign uOn, and correction to the vertex 
function equals 



r^if(n, n') = E ^('^ ~ n")Sn'+n"~nSn"a'n'+n"-nan"- (19) 



■n" 



One can find that 



r^^^(?T,, n') = 1 — A-j-fl7"(7r — 2| arctan arctan 



4 W Qph Qph 



and so does not depend from q. However, for g ~ which give the main contribution to 
the self-energy it presents the correct order of the correction to the vertex function (compare 
with Eq.(|^)) and the self-energy. From Eq. (2) one has for u/Qph ^ 1, Qph/W <^ 1 

This is similar to Eq. (|1^, but has a different numerical coefficient, due to neglecting 



the g-dependence in the vertex function ([T9|) 



V. C. GRIMALDI, L. PIETRONERO, AND S. STRASSLER'S APPROXIMATION 

In the paperlii the authors started with the same equation as in Migdal method, but made 
a series of approximate assumptions (for example, expanded in small g's to take integrals) 
and gave the following estimate for the vertex function: 

TQpg{uj, Q, f2, W) = -^[[arctanm — arctan ^ ^ — ] arctan — 

-[Q-m arctan ^] ^ ^ 7,,, ' ] • (20) 

Here they set one of the external electronic frequencies to zero. In this case this result 
is similar to that in the Migdal's approximation (|^. However, to calculate S one needs 
dependences on both external frequencies. 

There are some indications (see, e.g.0) that for a small hole doping strong Coulomb 
correlations renormalize the EPI, giving rise to the strong forward (small-g) scattering peak, 
while the backward scattering is strongly suppressed. With this idea in mind, in Ref.0 
the electron-phonon coupling constant was assumed to have a cut-off in a momentum space 
Idpkl"^ = g'^{2kF/qc)'^9{qc — \p — k\). An approximation for g was used@ and then 
Qc = Qc' 2pF was set to q^ = 2pp giving the following result^! for Z: 



'CGPS 



1 + X--X' 



4 W Q 
For = 1 it becomes 
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ZcGPS ^ 1 + A - O.SA^^, 

which gives the correct order of the correction to the self-energy , but underestimates it 
(compare with Eq.(|l7D). Probably this is due to the expansion in small g's. 

VI. TAKADA'S GISC METHOD 

According to Takada's gauge-invariant self-consistent (GISC) method, the vertex func- 
tion can be chosen as a functional of the self-energyi Vt = Vt^t]- This choice is based on 
the Ward identityi (^ valid for all ojyS and g's. In Ref-i Takada proposes neglecting the 
vector term in the Ward identity and choosing!! 



2iuJn 2{iUn — iujv) 

which corresponds to the estimates of RefsBi in this limit. If this method worked it would 
allow one to calculate also superconducting diagrams without complicated integrations over 
momenta. 

The set of equations (1,2) is suggested to be solved by iterations. At the first step 
G = r = T^^'^ = 1 give St = T.^t^ = S^^) from Eq.(2). At the second step G = G[I:^t^] 
; and the correction F^-* for ilph/W <^ 1 gives0 



Sg)(.c.) = -^X'n,, { + '''' . (21) 



TT^ ^ph I, ,1 ^ O 



For small frequencies it results in 



Zt = 1 + X + X\'^ + ^) ^ 1 + a + 0.9A^ 
lo 4 

which has the incorrect sign and order of magnitude of the vertex correction (compare with 

(O))- 

(2) 

During this calculation, the correction to the vertex Fj, is of the order of A for any q, 
which contradicts the Migdal's theorem! and can be valid only for g ~ (see e.g. Ref.i). 
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So in Takada's method the correction does not have the small parameter Vtph/W . This 
can be explained by the fact that the region Vf|^ ^ |^| where F ~ A gives in fact a small 
contribution to S. If, on the other hand, Vplg] ^ |z/|, the vertex V is of the order of XVLph/W , 
which gives ~ XnlJW (see (|13|)). 

To explain this result, one can find the vector vertex function F in the lowest approxi- 
mation and show that it really cannot be neglected. The diagram in Fig. 1 corresponds to 
r(2) jf p(i) jg substituted by F^^^ and F(^)(zu;„, iojn — iuJu, ^ - by T^'^^ito, ito — iu, —q). Using 
electron-hole symmetry, in our model T{iuj^ icu — iu,k,k — q) = T^^^ + F*^^-', where F*^^-* ~ ^, 

1 /•! ^ nQr] 



qf^'^\iuj,iuj — iu^q) = X- [ drj 

4 J-i 



—im' + Qr] 
1 + im + I — ifn I ^ j^ui — im' 



x{ln ^ + ln- 

^ l-im + — 1 



+ ln 



im + — 1 + im 1 — im + im' 

mo 

1 — im + im' — Qri + — 
+ ln — (22) 

1 + im-im' + Qri + 

Neglecting the Q-dependence in the last term of (p2D one can get 

gf = —iXflph(l — arctan ^^)(arctan arctan ). (23) 

We can find regions where scalar or vector terms in the Ward identity dominate. 

If one uses at ^Ip^/W <^ 1 the expression (0) for F^^^ one can find that the terms i/F*^^) 
and qT^'^^ from (|1]) are of the same order if u/qVp = 0.43. If, however, vjqVp < 0.43, the 
vector term dominates and cannot be neglected. 

VII. F. COSENZA, L. DE CESARE, AND M. FUSCO GIRARD'S 
IMPROVEMENT OF THE GISC METHOD 

Taking into account a criticism of the GISC method0, F. Cosenza, L. De Cesare, and 
M. Fusco Girardii tried to improve the GISC method and did not neglect the vector term in 
the Ward identity. They employed the simplest choice of the solution to the Ward identityii 



{iujn' — iuJny — \a{k, k')\'^\k — k 
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where a{k,k') = {q/2 + k)/m ~ Vp- With S(A;,zti;„) taken from Eq.(|^) one gets for |g| 
\k — k'\ ^ Pf 

T§l)G{k,iuJn,k\iuJn') ~ K-^y 

which is different from the Migdal's estimationi r*-^-* ~ (see (|^). 

Substituting approximation into the equation for the self-energy (§) one can get 



sgL(^^-0)-0.25.AM^f(ln^)^ 



which is different from (|l^) and gives 



ZcDG ~l + A-0.25A^(-^)^(ln , . 

This can be explained by the fact that the Ward identity is a one equation for two 
functions (the scalar and vector vertex functions) and thus allows for multiple solutions^. 
So there is no particular preference in using (^) instead of any other solution. 

VIII. DISCUSSION 

So in the paper we considered consistently the contribution of nonadiabatic effects to 
the electron self-energy in the normal state at T = in second orders of the EPI constant 
A. Several methods of taking into account nonadiabatic corrections were compared which 
can be summarized in the Table {pfVf = W). 

The results following from the approximations of V. N. Kostur and B. Mitrovidni, E. 
Cappelluti, C. Grimaldi, L. Pietronero, and S. StrassleJl^ give correct estimations of the 
order of magnitude and sign of the nonadiabatic corrections to the self-energy. The analytical 
results were compared with numerical calculations. 

It is shown that the GISC methodic and its generalization^ give for the self-energy 
overestimated and underestimated results respectively. This is connected with the fact 
that the leading contribution to S comes from the region qvF ^ v where F ~ XVlph/W . 
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This means that in the Ward identity one may not neglect the vector term as proposed by 
Takada. In the framework of the standard (Migdal's) approach, in the lowest order, there 
is a calculable parameter v/qvp ( the ratio of the frequency of an incoming phonon to 
its momentum) for which the contributions of the vector and scalar terms become of the 
same order. They are equal when v/qvp = 0.43. So for the most metals, where the main 
contribution to the self energy comes from the region v/qvp ^ 1, the GISC method may 
not be applied. One should also be careful when choosing a solution of the Ward identity, 
so as not to underestimate the scalar term, like in Ref.ii. It seems that GISC method can 
work in systems with long range interaction, for example, in doped semiconductors. 

The difference in the renormalization factors Z in the Table are due with the different 
approximations used in the calculations of the momentum and frequency dependence of 
vertex function. We can illustrate this in Fig. 3 where we represent the results for F*^^) 
in the Takada's (g-independent) approximation!, the Migdal approach for different q (Eqs. 
(4,27) of the present paper), and the g- independent Kostur, Mitrovic vertex function. We 
see that if the former overestimates the vertex corrections, the latter underestimates them. 

For a stable system the condition ImE(ict;„) < for cu^ > must be satisfiedS (it 
is equivalent to Z{iujn) > !)• However, as one can see from the Table, the nonadiabatic 
corrections reduce the renormalization function Z, favoring the tendency towards instability. 
The critical Xcrit at which it takes place (see Eqs. (10-16)) is given by 

arctan m — arctan 

\ _ l/mo+l 
■^crit 



mody{m, mo) + dr{m, mo)) 

for m — > 0, where d^{m,mo) and dr{m,mo) were defined above. The dependence of A on 
the parameter Qph/W is shown in Fig.4. Thus, for a stable system there are constraints on 
A for a fixed Migdal parameteJ^. This can explain, for example, the existence of PbBi with 
A ~ 3. At the same time, in some papers, e.g.0, it is stated that A ~ 1 independently from 
the Migdal parameter. As one sees from Fig. 4, it may be the case only for Qph ~ W. 

This possible violation of the analytical properties of the one-particle Green's function 
could result in corresponding nonanaliticity of the two-particle Green's function, and thus 
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in a charge response function. This can lead to a charge instabihty. 

It follows from our analysis that different approximations to the vertex function give 
in the normal state quantitatively (and sometimes even qualitatively!) different estimates 
for the self-energy. This makes the conclusions of these theories doubtful, as the same 
approaches were used there to calculate in the superconducting state. This means that 
only direct Migdal-type calculations can give an answer to the question: Do nonadiabatic 
effects enhance critical temperature of the superconducting transition?ii 
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^(2) 


A/TnTnpl ( mTmPTir'pl 1 wifn T^v / tttuti r^,n (4i 

IVJ-lgjU-CXi I llLLllid H^CXi 1 ^ VVluli J. 11 Will J_JU.l^) 


—9 TA^O 1, IW 


A/TlOTiPll ( fl Tl fl 1 A/'fl P^l 1 1 ■^X/'l"^n T^V / TTOTTI riH l7l 
IVl-lgjU-dl I cXllcXl^ IjlL/CLl ( ^ Wlull J- 11 vjlll J_Jl^. I I J 


—2 TA^O . IW 


r!T,X Rpf@ 

v^yuyv^ J. VjCI .1 — 1 


—1 '^SA^O u IW 

-L.OO/\ ^^phf 


KM Refii 


-1 73A^fi 


CGPS, Ref.0 




Takada's CISC, Ref.i 


+0.9A2 


CDG, Ref.i 


-0.25A2((]p;,/iy)2(lnTy/np;,)2 
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Figure Captions 

Fig.l. Equation for the vertex function T^"^^ in Migdal's method. The outgoing hues are 
shown for clarity and are not included in the definition of F. 

Fig. 2. a) Plot of F(2)(cj = 0,iu) (Eq.(4), dashed line) and Fg),^(cj = 0,iu) (Eq.(7), 
solid line) for A = l,nph/W = 0.1, u = 0, q/pp = 0; 0, 2; 0, 4; 0, 6; 0, 8. b) Corresponding 
frequency dependence of T,^\iLj) (dashed hne) and T.^^appr^'^K'i^) (solid line). 

Fig.3. F(2)(cj = 0, z/)in the Takada's approximation! (long dashed line); the Migdal 
approach for q/pp = 0, 0.5, 1 (solid lines, from top to bottom) ; and the Kostur, Mitrovicil 
vertex function ( short dashed line ). 

Fig. 4. Xcrit at which J]{iuj) = vs the parameter Qph/W. 
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